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The  problem  treated  1st  Find  the  optimum  pool  si.se  where  a pool 
of  vehicles  serves  the  needs  of  a group  of  people  and  equipaeuf  nay  hi 
rented  and/or  owned, 

The  planning  period  under  consideration  is  1 year,  which  if; 
divided  into  periods.,  Demands  cose  from  separate  distributions  for 
each  period,  Three  distinct  ways  of  supplying  equipment  are  of  fattrps  , 

1,  Owning  all  equipment  • 

2,  Renting  all  equipment « 

3o  Both  criming  and  renting  equipment. 

The  number  cf  owned  vehicles  is  held  constant  throughout  the  year, 
but  the  number  of  rented  vehicles  night  vary  frcea  period  to  period 
Once  an  optimum  schedule  is  determined,  changes  are  not  allowed  :.ii 
the  end  of  the  year  at  which  time  a new  optimum  schedule  may  be  de-set 
mined.  Each  pool  of  equipment  is  treated  as  a multichannel  queue  d 
system  with  parallel  channels— a common  waiting  line.  Demands  Do 
equipment  are  treated  as  arrivals  to  the  queue  and  waiting  time  r 
in  shortage  costs.  Waiting  time  io  treated  as  a function  of  the  nr  . r 
of  channels  „ 

Suppose  the  optimum  number  of  owned  vehicles  for  some  year  is 
smaller  than  the  optimum  for  the  preceding  year  a It  is  reasonable  oo 
assume  that  the  difference  will  not,  be  large  and  can  be  taken  care  of 
by  emitting  to  replace  seme  of  the  vehicles  already  scheduled  fc_  rep: 
meat.  Accordingly  we  begin  the  1-year  planning  period  at  the  time  of 
replacement  and  assume  vehicles  are  not  replaced  at  any  other  time,, 
Finally,  we  assume  setup  costs  are  incurred  only  at  the  time  vehicles 
are  rented  and  that  both  rental  and  ownership  costs  are  proportional 
to  the  number  of  vehicles  for  each  period. 


UITEOXWinXi 


ilpsKmt  Seeds 


Precise  tools  are  needed  to  determine  whether  organisations  hav 
too  much  equipment  or  too  little , Benefits  from  employing  equip  -oat 
should  be  balanced  by  the  costs-  of  providing  the  .equipment;  too  raich 
equipment  results  ta  excessive  ownership,  rental,  and  setup  costa 
whereas  too  little  eqiaipmeat  results  in  excessive  shortage  coats,, 

A method  la  needed  for  determining  the  best .balance „ 

The  traditional  ways  of  determining  how  much  equipment  is  needed 

ares 

—Arbitrary  assignment  by  published  criteria  (e.g*,  Assignee:  t 
standards  in  handbook) 

—Standards  of  utilization  (minimum  period  of  usage  or  minim  i 
mileage ) 

—Judgment  of  management  (intuition) 

The  assignment  of  vehicles  to  an  individual,  for  Ms  exclusive 
primary  use  Is  also  ofteu  traditional. 

In  this  paper  we  shall  discuss  a way  of  determining  haw  much  fleet 
equipment  should  be  owned  and/or  rented* 

Some  of  the  traditional  methods  of  measuring  utilization  are  as 

" 

follows: 

lo  Express  the  hours  of  use  as  a percentage  of  total  hours  a -ail:  blea 
2o  Count  those  day3  which  show  any  sort  of  use  and  express  he 
count  as  a percentage  of  total  available  days* 

Both  of  these  methods  have  shortcomings  because  they  portray  only 


"service;1'  i0e»,  they  show  only  when  a piece  of  equipment  was  buss1. 
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^e  equipment  was  t usyy 

Focus  on  Fleet «Typ<s  feaTyais 


In  this  study  we  focus  our  attention  on  equipment  which  can  serve 
the  needs  of  more  than  one  person  or  functional  group,  ¥e  call  such, 
a group  a pool,,  The  tf»  8,  Forest  Service  operates  large  numbers  oi 
such  pools  toro ughouro  the  United  States,  These  pools  will  provide  t'ue 
primary  examples  for  this  paper.  Some  examples  of  Forest  Services  pools 
are;  firefighting  aircraft  within  appropriate  regions;  sedans  mii  hick- 
up  trucks  for  those  headquarters  which  are  big  enough,  to  warrant  ny\‘\  •»{■>» 

' r " 

of  equipment,  as  pooling  is  defined,  The  reason  for  including  tvo 
equipment  «*ypes  xs  that  <#he  technique s of  analysis  are  similar  for 
studying  both  a .cleet  or  aircraft  and  a fleet  of  vehicles.  What  - ?'■ 
learned  from  studying  one  fleet  will  be  useful  in  studying  the  other. 

Statement  of  Problems 

Tae  tools  of  analysis  we  prepose  to  develop  are  to  be  used  ox-,  the 
following  type  of  problems  for  an  annual  planning  period: 

A*  Qptlmm  pool  size --all  equipment  owned . --Assuming  th  * set 
of  equipment  for  each  operating  pool  consists  only  of  owned  l. 

what  is  the  optimum  number  of  pieces  of  equipment  for  each  pool? 

*1  / 

Optimum  pool  sise— -all  equipment  rented.^/  »~When  dtmnd 
seasonally , if  the  number  of  pieces  rented  is  allowed  to  vary  doxi.  • 
the  year,  what  is  the  optimum  policy  for  scheduling  rentals? 


l/ 

- Included  within  the  definition  of  rented  vehicles  are  those  old 
vehicles  designated  for  replacement  which  are  held  over?  during  osa: 
periods  even  enough  replacements  have  already  been  received. 


~2 
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We  begin  by  dividing  the  year,  into  a finite  number  of  nerxoa  . , 
Assume  the  parameters  come  from  separate  distributions  for  each  ceric:* 
These  distributions  are  not  necessarily  similar.  We  have  chosen*  for 
simplicity,  a sample  of  a steady  state,  Poisson  Arrival- -Iferpor e rt-iul 
Service  System. 

We  obtain  shortage  costs  by  considering  each  pool  as  c murti-- 
ehannel  (parallel)  queueing  system  with  a common  queue,  each  v.ehi 
corresponding  to  a single  channel.  Demands  for  equipment  ere  treat el 
as  arrivals  to  the  system*  Each  trip  corresponds  to  a service  com.nl r 
•For  simplicity^  we  assume  the  vehicles  (channels)  have  the  sa 
mean  service  rate,  and  that  queue  discipline  is  first  come  first  i »■ 
There  are  two  distinct  but  equivalent  ways  of  arguing  the  queueing 
process— the  "limited -length  queue"  and  the  "servicing  of  machine s," 

r 

Both  lead  to  the  same  set  of  steady-state  equations.  The  former 
involves  the  idea  of  a "truncated1  arrival  process,  while  the  latte, 
does  not  and  therefore  is  perhaps  easier  to  relate  to  the  actual 
physical  situation.  Feller  (1957)  describes  a model  of  V machine 
serviced  by  "r"  repairmen  which  is  analogous  to  our  queueing  proc  t , 
we  label  the  machines  as  customer's,  the  servicemen  as  vehicles,  and 
a machine  breakdown  calling  for  service  becomes  a customer  demand  for 
vehicle.  Feller  gives  the  detailed  arguements  starting  with  first 
principles  which  we  will  not  repeat;  we  merely  restate  the  steady- 
state  equations  as  a beginning  of  the  detailed  queueing  discussion- 


-4- 


dnaly  is  qVVthV  queueing  .w&Zi  Obtain  equations  for 

the  expected  value  of;  the  wait  in  queue  for  ea eh  period*  'Sae  aathemati 


model  uses  a general  expression#  (X.  + X)  for  the  sasajat wait  in  oueue 

qi  1 

. 

for  the  ith  period  without  regard  for  Its  explicit  foray  therefore. 
the  probability  laws  governing  arrivals  said  services  need  not  be  the 
same  fox-  each  of  the  "m"  periods . 

Description  of  Costs 


Vehicles  must  be  modified  to  be  suitable  for  Forest  Service  use, 

' 

Ordinary  vehicles  are  not  suited  to  emergency  use  in  firefighting  a k 
for  operating  on  unsaved  roads.  Such  modifications  account  for  fcl. 


bulk  of  setup  cost  "A"  (dollars  ner  unit)  for  rented  equipment;  however 
the  omied  equipment  setup  cost  is  amortized  over  the  lifespan  and 

becomes  tsart  of  "B"  the  overall  ownership  cost  ner  unit.  Accord:'  ig  y 

- 

we  assume  setup  cost  is  incurred  only  for  rentals  and  only  when  they 
are  first  acquired,  not  when  they  are  returned  to  the  rental  agency, 
Otmership  cost  "B"  (dollars  per  unit)  is  computed  as  follows s Tr e 
total  costs  over  the  life  of  a vehicle  including  initial  Investment, 
setup,  interest  on  investment,  operations;  and  maintenance  are  divided 
by  the  total  number  of  periods  in  its  life  to  arrive  at  a series  of 
equal  payments  amounting  to  "B"  per  period.  The  purpose  of  this  - -er- 
simplification  is  to  put  rental  and  ownership  costs  on  the  same  basis  , 
By  so  doing,  we  overestimate  the  variable  costs  during  slack  periods 
and  underestimate  them  during  busy  periods.  Variable  costs  which  are 
approximately  the  same  for  rented  and  owned  equipment  should  be  c 1 Vi -a 
from  calculations  of  cost  coefficients  ’B"  and  "d^."  This  value  may  be 
recomputed  each  year  to  take  account  of  changes  in  the  size  of  the 
owned  fleet  as  well  as  changes  in  the  cost  items  mentioned  above. 
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Rental  ooot-  Md!*  por  unit-  includes  rental  agency 


end  costs  of 


operating  and  maintaining  the  equipment,  It  is  proportional  to  the 

number  of  vehicles  rented  per  period. 

Shortage  cost  {dollars  per  unit  time)  "K"  converts  risen  waiting 
time  fl(li)  into  shortage  cost.  Waiting  tine  means  time  spent  by  users 
who  postpone  their  trips  while  equipment  is  busy,  The  most  serious 
difficulty  in  matching  actual  characteristics  of  real-life  operations 
is  in  estimating  shortage  costs.  How  does  one  answer  the  question1 
What  is  the  cost  of  an  hour's  wait?  Could  the  waiting  customer  work  at 
something  else  meanwhile  or  was  the  time  completely  lost?  If  the  tri 
was  for  an  emergency  (such  as  firefighting),  the  cost  of  waiting  is 
much  higher.  One  would  expect  that  the  cost  of  waiting  varies  seasonally, 
being  higher  during  the  busy  summer  season.  Passenger  -carrying  vehicles 
are  employed  primarily  to  save  people's  time;  thus,  we  might  relate 
shortage  cost  to  the  wages  of  those  who  are  waiting. 

Prom  year  to  year,  the  characteristic  "mix"  of  the  fleet  may  change, 
necessitating  an  adjustment  of  the  cost  coefficients. 


MATHEMATICAL  FOBMJLATIQSf 


The  Objective  Function 

For  the  ith  period  of  the  overall  annual  planning  period, 

(i  55  1,  2,  •»«,  m)  let* 

d^  « Rental  cost  per  vehicle  during  period  ”i.” 

Si  « Humber  of  vehicles  rented  during  period  "i . " 

Y » Humber  of  vehicles  owned  during  the  year  (constant  over  all 
periods). 

= Humber  of  vehicles  available  (both  rented  and  owned)  daring 
period  "i." 


Mi  " Xi  + * 


W : 4 ) 

If  .£ 


Sxpeeteo  fciac  during  period  *’i’ 

are  available). 


■ *h?m  fc. 

i 


per  unit  time  of  wait  during  period  i. 
For  each  of  the  periods  let; 


2£r.?Xt 


A = Setup  cost  per  rented  vehicle. 

B «*  Ownership  cost  per  owned  vehicle  per  period. 
Then  for  the  ith  period; 


diXi  ~ Eental  cost  for  renting  X vehicles. 

BY  = Ofmership  cost  for  awning  Y vehicles. 

KiTV  Mi  = Shorta2e  cost  for  having  M vehicles  available 
(Mi  - X.  + Y). 

A maxCOjX^.^^)  = Setup  cost  for  renting  additional  vehicles 
XQ  » Initial  level  of  the  rented  fleet  at  time  zero. 


Total  cost  for  *'ith"  period  is; 

A ^ ) + d^X^  4-  K^F^ (M^ ) + BY 

For  "m*  ueriods  we  have  total  costs 
m 

Tm  ”1  [A  + aixi  + KjS  (M  )]  + 

ial 


BY 


For  the  case  of  only  rented  equipment,  we  set  Y - 0 and  M « X 


1 i 


giving; 


m 

V""*>  .£  [A  * KlVMi 

i=l 
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For  th.  ch--. 
giving*. 


' . only  j-ffls c.  •*vili»xsr».t  e.  s;et  all  X,  v-  0 •.  -ad  U,  «*  T 


El 


T 


m 


(own)  « ) [K.ft  .(X)]+  :mBT 
L-i  x qi'  ■ 


1=1 


IV/ 


Mean  Walt  in  Queue 
For  the  ith  period  (i  = 1,  2,  . „«f  ®)  let: 

?•».  Mean  rate  of  demands  per  unit  time, 

U,  » Mean  length  of  time  per  trip  (mean  service  time). 


N<  = Maximum  number  of  people  who  demand  equipment. 
f'L  * Humber  of  pieces  of  equipment  (channels ) avail 
W . = Mean  length  of  time  waiting  in  line. 


qi 

L ~ Mean  number  of  people  in  the  system  both  waiting  in  line  and 
occupying  vehicles. 

E = Mean  number  of  people  waiting  in  line  only. 

Pa  a Prob. (the  number  of  people  in  the  system  is  V). 

P = ~ 

Subscripts  ("i" ) are  now  dropped  since  the  steady-state  equations 
are  independent  of  the  period  of  the  year. 

Assuming  Poisson  Arrivals  and  identical  Exponential  Services  times 
we  have  the  following  solutions  to  the  steady-state  equations: 


p = p- 
n 0 


pn  ( ^ ) ....  ( 0 < n < M) 


n IV 


n; 


n o 


nJ  M,  Min 


(Min/N) 


m 


■8- 


L a Expected  value  of  a,  the  number  in  the  entire  gyateeu 


CO 


We  need  a convenient  expression  for  W the  mean  wait  in  over- 

q * * 

Rather  than  argue  from  first  principles,  which  is  tedious,  and  be 
it  is  convenient  to  use  later  in  a discussion  of  discrete-convexity  of 
the  objective  function,  we  will  use  a shortcut,  namely,  that  the  mean 
wait  in  queue  is  a linear  function  of  the  mean  number  waiting  in  queue 
Justification  for  this  may  be  found  in  Little  (i960).  Ee  shows  that 
if  the  three  means  L,  X,  W,  are  finite  with  corresponding  stochastic 
processes  strictly  stationary,  and,  if  the  arrival  process  is  metricc.il 
transitive  with  nonzero  mean,  then  L = }£, 

Where: 

L = Man  number  in  the  system. 

W » Mean  time  spent  by  a unit  in  the  system. 

- Mean  time  between  two  consecutive  arrivals  to  the  system. 

Following  Little  (i960),  we  define  the  system  as  the  queue  itself 
and  let  L ^ refer  to  the  mean  number  and  wait  in  queue  and  obtain 

<1  q 

£,  * X&  . Our  queueing  system  has  finite  means,  is  strictly  stationary 
and,  as  Doob  (1953  p.  ^60)  points  out,  is  metrically  transitive  because 
the  random  variables  representing  the  lengths  of  time  between  arrivals 
in  a Poisson  process  are  mutually  Independent  and  from  a common  dis- 
tribution. 


* 


OUCttt:  ill 
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L/l  as  the  ejected  wait  in  the  system  we  hares 

E [delay  in  the  system]  ~ I [delay  in  queue]  + E [delay  in  service] 

or: 


X 


$ 


1 


'8J 


Using  equations  (j)  and  (8),  vs  arrive  at  a convenient,  closed-fora 
expression  for  the  mean  trait  in  queue: 


-1  p0 
W * 4*  ‘~r~ 

q p X 


M 


S 


n*l 

y 


C.  \ . .iM 

Hi  if  3 
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It  may  happen  that  the  values  for  mean  trait  in  queue  obtained 
from  the  finite  queue  model  may  he  closely  approximated  by  the  infinite 
queue  model.  Since  the  infinite  queue  model  equations  are  more  tractable, 
it  might  he  used  instead,  Therefore,  ve  list  the  corresponding  equation 


with  a 


Bo* 


superscript  to  denote  that  Infinite  queue  is  allowed. 


M 


r 


q 


-i  po  p 
+ — — 


“ ’ « »t%)2 
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Maximum  queue  length  equals  (K  -M ),  the  difference  between  U the 

maximum  number  of  neople  who  demand  equipment  and  M . the  number  of 

. i 

channels  (vehicles).  When  the  number  of  vehicles  is  large  enough 
to  allow  one  for  each  user,  i.e.,  - Si,  then  no  queue  can  exist  and 

mean  wait  is  zero  (which  suggests  why  we  assume  a finite -length  queue 
model).  N,  may  be  larger  for  some  periods  than  others  and  an  optimum 
value  for  Y may  be  larger  than  some  3>h , smaller  than  others.  However, 
an  optimum  value  for  3L  may  not  exceed  for  to  do  so  is  bo  have  a 
surplus  of  vehicles.  Thus,  the  rented  equipment  (X, ) may  vary  in  number 
but  is  bounded  as  follows:  X.  4 N.  (i  = 1,  2,  m). 
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equations  (5),  (6),  and  (7)  in  which  we  shall  use  the  property  that, 
our  objective  functions  are  "discrete-convex; " ;Ue.,  it  has  non -nega- 
tive second-order-  forward  finite  differences* 


We  investigate  the  property  of  discrete-convexity  for  both  finite 

and  infinite  queue  models  in  the  appendix,  in  detail.  For  the  rcj^.dt, 

we  ore  concerned  only  with  the  property  that  for  both  models,  the  naan 

trait  in  queue  is  a discrete-convex  function  of  the  number  of  channels, 

l.e.,  (iL  ) =s  Y)  is  discrete -convex „ 

Next,  we  examine,  the  individual  terms  of  the  objective  function 

and  show  that  each  is  discrete-con vex.  Both  C.X.  and  eBT  are  linear 

i 1 

and  therefore  discrete -convex.  That  Max(0,X^-X^^)  is  discrete- convex 
is  not  so  obvious  that  a graph  would  not  help  to  see  this 3 


suxCOjX^ 


& 
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We  could  join  discrete  points  by  straight  lines  to  form  a continous 
convex  function  from  each  of  our  • integer -valued  functions.  Since  the 
positive  sum  of  convex  functions  is  a convex  function,  we  have  shown 
that  if  our  objective  function  "f"  were  made  into  a continuous  function, 
it  would  be  a convex  function  (with  the  desired  property  of  having  non- 
negative, second-order,  forward  finite  differences).  Finally,  restrict- 
ing X.;  and  Y to  the  integers  doesn't  change  this  non-negativity  and  we 
have  shown  that  our  objective  functions  are  discrete-convex. 


In  the  inventory  problems  treated  in  the  literature,  the  u*.nt.lty 

of  a commodity  13  usually  treated  as  a continuous  variable;  ben.ee,  /or 

convex  objective  functions,  it  follows  'that  a local  minimum  i also 
global  minimum.  This  is  also  true  for  the  multiple-period  model 
function  of  several  variables.  However,  in  our  case,  expected  wait 
in  queue  is  defined  only  on  integers.  Elsewhere,  we  shewed  our  cbje-v.  ive 
function  to  be  "discrete -convex”  (in  the  sense  that  second- order, 
forward  finite  differences  with  respect  to  each  variable  are  non-naga- 
tive ) o But,  a discrete-convex  function  of  several  variables  car.  Lave 
a local  minimum  which  is  not  a global  minimum.  We  will  slier;?  an  example 
in  the  appendix. 

Furthermore,  we  have  an  additional  complication  due  to  the  discrete- 
ness of  the  variables.  Suppose  we  could  redefine  our  c ;j  active  fu&o  ,icn 
as  one  of  continuous  variables  with  the  same  functional  values  when  ver 
the  variables  are  integers,  (i.e.,  define  a smooth  surface  through  the 
vector  points  in  m-space).  Then,  assuming  a convex  function,  we  still 
cannot  guarantee  that  having  found  its  absolute  minimum,  we  can  roar.  1 
off  the  continuous  variables  to  integers  and  still  have  the  s yluti.br 
to  the  corresponding  problem  which  was  restricted  to  the  integers. 

The  example  in  the  appendix  illustrates  this  idea  as  well  as  the 
property  that  a global  minimum  need  not  be  coincident  with  a local 
minimum,  for  a discrete-convex  function.  For  this  example,  the 
author  is  indebted  to  his  colleague-  Gherman  J.  O'Ueill. 

Further  analysis  is  necessary  to  check  if  the  equation  for  oar 
objective  function  exhibits  the  "trough-like"  characteristic  of  this 
example  and  if  so,  whether  in  any  real  case  it  might  not  give  a good 
approximation  to  the  integer-valued  problem. 


*. 


For  'both  convex  and  discrete-convex  functions,  the  &5.ni,  urn  jiay 

lis  on  a boundary  of  the  region  on,  which  the  functions  are  defined 
rather  than  at  an  extremum,  /my  solution  procedure  East  take  this 
into  account. 

However,  the  situation  is  not  hopeless.  For  the  case  of  owning 
all  equipment,  we  have  only  one  variable,  "Y**  and,  for  discrete -con- 
vex functions  of  one  variable,  the  foregoing  difficulties  don't  exist 
because  a local  minimum  is  also  a global  minimum.  Per  th>.  cases  of 
(a)  rent  on3_y,  and  (b)  rent  and/or  own,  we  shall  introdw.ee  dynamic 
programming  to  overcome  the  foregoing  difficulties,  check! <g  boundary 
points  as  well  as  achieving  an  absolute  minimum. 

ALGQRFXHM3 

Dynamic.  Progx^aEming  Solution  rf:  Rent  and/or  O. Problem 

Rather  than  discuss  a separate  algorithm  for  the  case  of  rent 
only,  we  shall  discuss  the  more  general  situation  of  rmt  and/or  buy, 
treating  rent  onlv  is  a special  ease  when  Y = 0.  Using  the  recursion 
methods  of  dynamic  programming  seems  appropriate,  because  the  variab3.es 
X^,  ih,  and  n are  integers  (See  the  preceding  discussion  cu  convex 
and  discrete -convex  functions  of  severa3.  variables..  Furt j. t rraore 
these  recursion  methods  lead  to  an  absolute  minimum  (Belluai  find  Dreyfus^ 

1962)o 

The  key  to  a dynamic  programming  formulation  of  our  problvn  lies 
in  developing  a recursion  relationship  between  two  adjacent  periods, 

"i"  and  "i-l."  This  is  called  a "functional  equation."  Ue  begin  with 
a one-period  model,  then  with  a two-period  model,  ano  generalise  t \ 
results  to  an  m-period  model. 
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Gas-Period  Model 

As  before,  X is  the  number  of  owned  vehicles*  For  the  moment, 
consider  Y as  a fixed  number.  Since  the  maximum  number  of  vehicles  we 
could  supply  is  bounded  by  the  maximum  number  of  people  the  vehicles 
serve,  we  have  the  restrictions: 

0 £ Y£N 


0£X.<  EL£  V 
i i max 


max 

(i  ™ j-j  2,  m) 

where:  XI  = xsaximics  value  any  N.  may  take  on  for  any  of  the  V 

periods,  i.e.,  = max^j-  (i  « 1,  2,  in).  Let:  X()  = initial 

size  of  rented  fleet  at  time  zero,  (X„  * 0,  1,  2,  .....  N ).  The 
total  cost  for  a one -period  model  is: 

Xq,Y)  a + A maxCOjX^-X^)  + Y)  «*  31 

We  have  designated  I as  a function  of  Xq  and  Y as  well  as  X^ 
because,  although  we  treat  Xq  and  Y as  fixed  for  now,  later  we  will 
let  them  vary  and  the  dependence  of  our  answer  on  X^  and  Y should  be 
kept  in  mind.  Minimizing  T(X^,X0,Y)  with  respect  to  yields  a new 
expression  C^(Xg,Y)  which  we  designate  as  a function  of  X„  and  Y with 
the  understanding  that  for  the  time  being  X^  is  the  only  variable: 

C-^X^Y)  = MnjTCXpX^Y)] 


where: 


Thus: 


CX(X0,Y)  = min [d^  ♦ A saxCO,^-^)  + l)  * BY] 


(n) 


(32) 


(13) 


JA- 


» 


Setting  Y » 0 gives  us  the  expression  for  the  case  of  rent  only, 
namely: 

Cl(V0)  . ^ ♦ A -(0,  VV  ♦ hW.  1 

(hr 

M-Feriod  Model 

Consider  the  case  of  a two-period  model 0 Renumber  the  periods.,  such 
that  the  first  period  in  time  follows  the  second  period  along  the  time 
axiSo  A graph  of  a sample  two -period  situation  may  help  to  see  ti 
relationships  : 


Period  2 begins  with  an  initial  rented  fleet  size  of  XQ;  however,  the 

initial  rented  fleet  size  for  period  1 is  Xg«  Beginning  at  time  zero 

A 

with  Xq,  again  holding  Y fixed  and  rent  Xg  = Xg  vehicles  for  period  2, 
we  want  to  minimize  the  cost  for  the  remaining  per' od  1,  whit h by 
definition  is  C^($g,Y)j  so  our  cost  for  renting  Xg  - X„,  would  be 
defined  as: 

Tg(Xg,X0,Y)  - dgXg  + A mx(0,X2-X0;  + Kg^g(Xg+  Y)  + C^X^q) 
for  a two-period  process  ■> 
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i-.  opi/iiiku.  choice  for  Xg  is  one  whxcit  idni&igcs  this  function! 
accordingly,  letting  again  become  a variable  ve  can  defines 


C2(X0?Y)  - ain  T2(X,vX0,Y) 


X, 


C2CX0,T)  » ^ittldgXg  + A ffiaxCo,X2-X0)  + K£f  ,;.(Xp+l)  4-  BY  4-  C (X^Y)] 


X 


•where: 


2}  * {X2 


V * fXg  Od^Bfg-Y 


It  follows  iaasdiately  from  the  preceding  two-period  discussion 
bnat  t,he  general  functional  equation  for  six  n?.—neriod.  model  is: 

Ca(X0’Y)  * A mas<°>VX0  + EAm(VY)  + + Vl(Xi’?J 

IV 

where:  jx^J-  * -jx^J  0 H^xj- 

(i  1,  2,  0000^  m3 

Bant  Only 

As  before,  setting  X - 0 we  have  a Multiperiod  Bant  Only  Model, 
namely: 

0^,0)  re  MnJdJC^  4*  A a>ax(0,Xm»X, 


XL 


m 


m n 


» an  0 


E # (X  ) + C ,(X  ,0}] 
m cp  m m-1  nr 


For  a fixed  X,  we  have  shown  how  to  find  Repeating  the 

procedure  for  each  value  assigned  to  X in  which  (0  € Y $ IT  )» 

••  0 % max'* 


(IS) 


(ill 


fi 


(19 


i2Cl 


l6“ 


■A 

- 


That  C (XqjY)  is  discrete-convex  is  not  obvious . To  them  this,  we 

argue  inductively  starting  with 

C-jfoY)  » min.  TfoX^Y) 

Wr 

G1(X#Y)  m minCd^  + A njax(0,X1-X)  * Kjff  ^Jt+Y)  + BY] 

IV 

We  have  already  discussed  the  discrete-convexity  of  T(X,XpY  for  any 
given  value  of  X^.  How,  we  show  that  T is  discrete -convex  in  X for 
fixed  (XpY),  and  is  discrete-convex  in  Y for  fixed  (X^X}. 


For  a given  X^  maxfO^-X)  « }X1“q  otherwise}"  11X13  is  graphed 


(13) 


below: 


MaxCO^-X) 


If  we  joined  the  discrete  points  by  straight  lines,  the  result 
would  be  a convex  function  with  non-negative,  second-order,  forward 
finite  differences;  and  restricting  that  convex  function  to  the 
integers  doesn't  affect  that  non-negativity.  Hence,  laaxCo,^  -X)  is 
discrete-convex.  We  have  already  discussed  the  discrete -convexity  of 
W where  = X^+Y,  so  that  for  any  fixed  X^,  ft (Y)  is  discrete- 
convex.  "BY"  is  linear  and  thus  discrete-convex.  Finally,  treating 
each  of  the  terms  in  (13)  as  a function  of  two  variables,  X and  Y,  so 
that  they  are  all  defined  on  the  same  domain,  we  may  add  together  dis- 
crete-convex functions  to  obtain  another  discrete-convex  function, 
namely,  C1(Xq,Y). 


We  turn  now  to  examine  the  discrete -convexity  of 


Gg(X0,Y)  - Brin  TCZg-^Y) 
X. 


fl  X 


where 


X. 


0*2^*ffg-Y 


Examining  the  last term  of  the  above  functional  equation,  for  fixed 

x2,y, 

C^X^Y)  « min  [drX(  * A raax(0,X- -Xg)  4 I^f  (X^  + Y)  4 BY 

M ‘ ~ 1 ; 

which  is  an  equation  of  the  same  form  as  C-,  (Xq,Y)  of  which  the  discrete- 
convexity  arguements  have  already  been  given,  ‘the  remaining  terms 
are  also  of  the  same  form  as  for  C3  (X^,  Y),  so  the  previous  arguements 
again  go  through , 

Finally,  assuming  that  ^ (Xq,Y)  is  discrete-convex  we  show  that 
Cn(X0,Y)  is  also  discrete -convex  by  again  referring  to  the  previous 
arguements  for  the  case  of  C^Cx^Y), 


lb; 


miafdgXg  4 A max{0,2g~X)  4 ^q,Pv  * t) 

Iv 

+ by  4 CjCXgjY) 


(IT) 


Computation  Procedure 

1.  Let  X be  the  number  of  rented  vehicles  at  the  beginning  of  the 
highest  numbered  period.  Compute  for  each  value  of  X *=  0,  1,  2,  . If 

Dicier 


each  value  of  X « 0,  1,  2,  , N,  and  Y = Xv"  a table  of: 

Y1(X,X1,Y<e>)  a [d1X1  4 K1Wa3(Xl  4 YW)  4 BY"V4  A maxCOjXj-X}]. 

From  this  table,  select  min  V1(X,X-,Yt‘' ) * C;l(a,Y^)  for  each  value  of 

|ci(X,Y^)| 


X giving  a table  of 


draw  to  evaluate  C^(Xgj?Y^ )/ 


values  from  which  we  may  subsequently 
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. S*JF 


* 


2*  Compute  for  each,  value  X * G,  1,-2,  . X , each  value 


of  X, 


0,  1,  2, 

AO 


mr 


and  again  fix  Y = , a table  of 


Vg(X,X2,YW  ) * [d2X2  4=  K2^2(X2  + Y4°)  + BTCU  4 A m3x(0;,Xg~X}  ] 

For  each  value  of  Xg  we  add  C (3^,  Y’^  ) from  step  1 giving  another  table 


r(0 


OKS 

c1(x2,y  " ) 4 v2(x,x2,y^3 

for  all  values  of  X from  which  for  each  X we  can  select  the  minimum 
value  giving  a table  of  C0(X,Y4*^)  from  which  we  may  subsequently  draw 


to  evaluate  CgCX^Y4^)* 


Simiarly  we  continue  until  we  have  finally  evaluated  Y*w)o 

Throughout  the  preceding  discussion,  we  kept  Y fixed  at  Y =>  Y '*  To 

. 

complete  the  solution  procedure  we  repeat  the  foregoing  for  all  values 
of  Y over  its  range,  (OjY&H^y)  in  which  * max^N^  For  each 
value  Y assumes,  compute  the  optimal  values  X^(Y),  Xg(Y),  X^Y) 
together  vlth  the  aeeodated  Ca(X0,Y)  and  from  such  a table  We  can 
select  the  final  answer.  Rather  than  store  the  results  for  each  value 


C«) 


chosen  for  Y,  we  can*,  compute  C^X^Y1'*)  and  C ) for  Y = Y 

t*> 

and  Y » Y rensectively;  compare  the  two,  discard  the  larger,  store 

. 

the  amallei*;  compute  Cm(XQ,Y4  ) for  Y * Y4*^  compare  this  with  the  pre- 
viously stored  value,  discard  the  smaller,  etc0 

So  far  we  have  only  briefly  discussed  how  to  evaluate  C,  (XQ  Y) 
for  any  given  i =1,  2,  m and  any  given  Y.  Since  Ci(X(),  Y is 

integer-valued,  methods  of  calculus,^  convex  programming  are  not 

applicable o A thorough  discussion  of  this  problem  lies  beyond  the 

- 

scope  of  this  paper.  (An  optimal  search  technique,  using  Fibonacci 
numbers  for  locating  the  coordinates  of  a maximum  for  unixaodal  functions 
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is  discussed  in  Bellman  and  Dreyfus  (1962,  v,  152)).  However,  for 

problems  of  the  size  we  expect  to  encounter  (12  periods, 

the  following  more  simple-minded  method  seems  adequate.  For  example, 

in  C-JX^Y)  » minCd^  * A max^^-X^  + K^f^CXj  + Y)  + BY]  we  have 

\ 

> 

a simple  linear  relationship  in  each  term  except  for  shortage  co3t 
^l^ql^l  * which  is  a function  of  ^.(X^  + X)  the  mean  wait  in  queue, 
For  convenience,  we  repeat  this  latter  expressions 


W 


ql 


(\  + Y) 


1^1 

n=l  n^+Y+1  MiI  wj 

This  is  a closed  form  and  can  he  computed  beforehand  for  appropriate 
values  of  Xi  + Y,  keeping  in  mind  the  ranges  0 4 (X^+Y ) g mexjV 
(i  = 1,  2,  . m).  The  remaining  terms  d^X^,  BY,  A raaxC^X^-Xg) 
are  easy  to  compute  and  store  so  that  finding,  a minimum  is  straight- 
forward,, 

Referring  to  equation  (2),  the  original  objective  function  of 
m + 1 variables  has  now  been  reduced  to  m discrete-convex  equations  of 

§ ’ ' *i  " l.y  ~ ; $ K 1 1 i T -• 

one  variable  (for  each  value  of  Y = 0,  1,  2,  II  ^ ),  Hence,  at 

each  stage  the  minimizing  procedure  guarantees  an  absolute  minimum 

t 

may  be  found. 

Alternate  Algorithm  for  Own  All  Equipment  Problem 
In  this  case  since  no  equipment  is  rented  we  set  all  = 0 
our  problem  becomes;  minimize 


Tm(Y) 


mBS  + 


IK!V(Y) 


i=l 


(9) 


(21) 
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where  Y is  an  integer.  Defining  Y*  as  the  value  of  Y which  minimizes 
T,  and  since  5?  (Y)  is  discrete-convex,  we  have  the  following  conditions 

Hi 

which  Y‘’:'  must  satisfy: 

A T<Y*>60 
A T(Y*+l)»  0 

where  A T(Y*)  is  the  first  forward  finite  difference  at  Y*. 

finding  Y*  satisfying  (22)  involves  a search  technique  of  which 
there  are  numerous  varieties.  Since  T(Y)  is  discrete -convex,  thi3 
property  can  probably  he  employed  to  find  a more  optimal  search  tech- 
nique than  the  one  following.  For  problems  of  the  size  we  expect  to 
encounter,  namely  where  YclOO,  what  follows  appears  to  be  a workable 
method. 

To  find  the  Y*  which  satisfies  these  equations,  assume  X*  is  unique, 
then  begin  with  Y = 1,  and  continue  increasing  Y by  one.  At  each  step, 
solve  for  the  finite  differences  indicated  by  equations  (22).  If  at 
the  onset  A T(Y  = l)>  0 then  optimality  has  already  been  reached 
and  no  further  steps  are  necessary. 

The  procedure  outlined  above  will  produce  a global  minimum  at  Y* 
since  for  discrete-convex  functions  of  one  variable,  a local  minimum 
has  the  same  functional  value  as  the  global  minimum.  If  Y*  is  not 
unique,  equations  (22)  would  read: 

A T(Y*) < 0 
A T(Y*+l)  « 0 

which  indicates  that  the  choice  between  Y*  and  (X*+l)  is  optional 
from  .a  purely  economic  viewpoint.  Hence,  a decision  is  based  on  other 
qualitative  grounds.  Accordingly,  we  have  shown  how  to  find  the  optimum 
pool  size  when  all  equipment  is  earned. 


CONCLUSIONS 


Referring  to  the  cost  of  wait,  for  the  Forest  Service,  these 

"values  of  postponement"  can  be  collected  as  part  of  the  overall  data 
collection  scheme  and  could  constitute  a wide  survey  which  in  Itself 
would  be  of  value  to  the  management  of  any  Region  wishing  to  analyse 
vehicle  usage  of  its  National  Forests  or  Ranger  Districts, 

Following  are  3cme  possibilities  for  obtaining  costs  of  waiting 
(postponement).  We  employ  aircraft  primarily  as  a firefighting  weapon 
and  a shortage  of  aircraft  results  in  increased  fire  spread j thus,  we 
might  relate  shortage  cost  to  the  effectiveness  of  aircraft  in  fighting 
fires  and  the  cost  per  unit  time  of  fire  spread. 

If  all  methods  of  estimating  shortage  cost  fail  to  produce 
reliable,  absolute  values,  we  can  still  use  these  estimates  in  studying 
renting  of  equipment  compared  with  owning  It  because  the  saortage 
costs  are  the  same  for  both.  We  can  also  let  waiting  cost  per  unit 
time  be  a variable  and  study  the  sensitivity  of  the  model  to  changes 
in  waiting  costs. 

Referring  to  "B"  and  . the  ownership  and  rental  costs,  a 

more  accurate  approximation  would  result  in  separating  each  of  these 
costs  into  two  parts,  fixed  costs  and  variable  costs.  Then,  variable 
costs  could  be  charged  on  the  mean  number  of  vehicles  in  use  as  well 
as  the  mean  service  time  (trip  length).  However,  for  this  study  we 
use  a fixed  coefficient  "B"  as  a first  approximation. 

Referring  to  the  assumption  of  steady-state,  most  motor  pools 
familar  to  the  author  are  characterized  by  having  trius  of  short 
duration,  with  a few  lasting  several  days.  If  one  lays  the  results 
of  each  day's  operations  end  to  end,  we  might  approximate  a realisa- 
tion of  a renewal  process,  an  example  of  which  is  the  Poisson  process. 


Although  the  dynamic  programming  functional  equations  are  formulated 
without  regard  for  properties  of  the  objective  function*  we  conjecture 
that  some  more  efficient  means  may  be  devised  to  perform  the  minimising 
operation  at  each  stage  of  the  calculations . The  reason  is  because 
the  objective  function  as  well  as  Cm(X^,Y)  are  discrete-convex  and  at 
each  stage  the  minimizing  operation  is  performed  on  a one  variable 
problem., 
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Proof  of  Discrete -Convexity  Properties  of  the  Mean  Wait  to  Queue 

TtrT — f.n. r u I,-,  w ip  »wr«iiiCT-iMi.i  ijfr  ■ in  i - itk  imiwnji)ini  u m i'-)i-i  -r  11  ng  i.rr  — i n-  rrr  -it  -t. n~  i- — i — rf-iy— tr- ‘v  — ir  — --r  — -rr  mnwiMw***1** 

Theorem:  For  the  Poisotei '•exponential,  finite -length  queue  system, 
the  mean  wait  in  queue  is  a discrete-convex  function  o ;’  the  number  of 
channels  available  during  any  period. 

We  first  show  the  following; 

Lemma;  Both  the  mean  number  in  the  system  and  the  mean  number  to  the 
queue  are  discrete  convex  functions  of  the  number  of  channels. 

Proof;  The  expression  for  £,  the  mean  number  in  the  system,  is  repeated 
for  convenience. 


For  convenience  in  what  follows  we  drop  the  subscripts  which  designate 
the  period.  The  second-order  difference  is  denoted  by 


£2L(M)  - p 


n=l  n*M4-3 


n; 


(MM)S  KH  f'} 


(Mp1  l 


Tliis  expression  can  be  simplified  by  letting 


and  noting  that 


Mj-2  M±1  Jl 

f (a)  + ^f(n)  s f{»2-2)  ~ f(M-i-l) 


n=l  n»0  n=0 

bo  that  £?L(M)  can  he  written: 


« f(M-M)  [1-2  -5*  j -J*  ft  M^*2) 

1=1  ‘"‘•{•± 


1-2  i (M>X 


3L 


f (Mti ) (K+i ) ! 

— "SI— “ 


(MU)1 


bj  1 


then: 

and: 

and: 

and: 


f(K*-i)fcO  for  (i  »1,  2 , I?) 


Mfl 

M 


>1  for  04  « I,  2.  . . . . ) 


M 


> 

1 


2||||j  for  (H  * 1,  2,  ) 


K * — rr  > -* — * it  * i,  a $ 

M (KM)(K*2)  (M+l) 

The  second-order  difference  is  non-negative  and  we  have  shown  that 

L(m)  is  discrete -convex.  Let  L be  the  mean  number  in  the  waiting  line, 

<1 

Frcra  equation  (8)  we  multiply  by  X,  giving: 

L ® X#  -J-  77 

<1  h 


or; 


L sl.p 

q 


Since  L is  discrete-convex  and  since  p is  a constant,  it  follows 


that  l 


is  discrete-convex. 


(27) 


(28) 

(29) 

(30) 

(31) 


(33) 


Again  from  Litte  (l96l)?  for  systems  that  behave  like  ours  w - Z /. 


q 


q/h 


so  that  ¥^(M)  is  convex,  since  it  is  a linear  function  of  which 
establishes  the  theorem. 


-25- 


Theorems  For  the  1 ’inf  Inite -length*1  queue  system,  tbs  mean  wait  in 
queue  (in  line)  is  a discrete -convex  function  of  the  number  of  channels 
(vehicles)  available  during  any  period. 

Proof; 

lets 

%'J  ~ Mean  wait  in  queue  (infinite  queue)* 

^ * -1.  % 


W 


q 


4 Mp-l 


Where! 


Q, 


Pa  M 

0 p 


"M  “ MS  „ p 

i.4«»rr 


M 


«° 


P0P 


P0P 


a 


fJL  t-f  #_ 


4fti 


U-g)  (Mu-X, ) 


Mi  14m(i  A }? 
Mi' 


(3't) 


(35) 

(36) 


Treating  this  as  a discrete  function  of  M,  we  show  that  the  second-order 
finite  differences  are  non-negative, 

A%°(M)  = f°(M5*2)  - 2V7°(M;-l)  -f  W°(M) 

» " a q 

M 2 

a2jsO/  , %p  p 

fil  W . fc  J “ «— — — n~r  |,M  , , . . 

* MiOft. S)(M+1)(1~j&-}  £<M-l-2)  u-X] 

M 2 

ft  V p 
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Since  each  of  the  three  terms  (without  regard  for  the  sign  before  the 
second  term)  is  non-negative,  it  is  enough  to  show  that  the  second 
term  is  dominated  by  the  third  terms 


(M-Kl)(l-^r)  [(Ms-1)  l ^ (l-§)  £M-X  3 


M' 


(39) 


To  see  this,  first  note: 


gfjr  * 1 ^For  p4lJ  M^2-) 

Then  compare  the  remaining  similar  terms  in  the  denoLu.ua tor 
First: 

4 r?  {for  p<l,  M*l) 


M*1  * M 


thus: 


Then; 


M+l 


>1 


_ £ 

M 


i*o) 


(l«£H)p  - X>Mp  - X (for  M^l) 

We  have  shown  that  W_.(M^),  the  mean  wait  in  queue,  is  a discrete- 
convex  function  of  = X^  4-  Y,  the  number  of  channels, 

A Discrete-Convex  Function  with  a Local  Minimum 
fsfine  a minimum  of  a discrete  function  as  follows : We  are  at 
a minimum  if  in  moving  to  any  adjacent  point  in  the  lattice,  ohe 
functional  value  increases  or  remains  the  same.  An  adjacent  point  means 
one  which  is  reached  by  making  a unit  change  in  any  or  all  of  the 
variables.  Consider  the  graph  showing  a lattice  of  a discrete-convex 
function  gCXpXg),  which  is  restricted  to  integer  arguments  of  a convex 
function  f (X-j,Xp),  defined  everywhere,  so  that  the  functional  value  of 
point  (X^  ,Y^)  i3  equal  to  the  perpendicular  distance  of  the  point 
(XpY^)  to  a diagonal  line  joining  two  nonadjaeent  points  of  the  lattice,: 
An  example  is  shown  belowj  the  dotted  lines  represent  contour 
lines  of  equal  functional  value,  f(XpXp)  * 1,  2,  3,  Only  a few 


(m 


m 


(^3) 


-£?~ 


contour  lines  .of  special  interest  hare  been  drawn » 53ie  surface  formed 


by  4f(>L  ,Sg)  is  a symmetric  "trough*  whose  bottom  fCX^I^)  ~ 0 lies 
along  the  diagonal  line.  Functional  rallies  are  shown  for  some  points 


indicated  by  square  brackets  on  the  lattice, 


Figure  3 


re  t± ne  both  functions  by  adding  the  plane  bfXpXg)  « s(3-X), 
s>0,  This  plane  intersects  the  X^Xg  plane  along  the  line  X1  = 3 and 
slopes  to  the  left  and  upwards  with  slope  «*  e>  0.  The  new  function, 

A + h'X^jXg)  « jCX^jXg)  for  small  we"  la  strictly  convex.  Restrict- 

ing it  to  integer  values  of  X^Xg  defines  a new  function,  k(X1,Xg)  which 
is  discrete-convex#  All  functional  values  of  kCX^Xg)  are  slightly 
different  from  the  corresponding  functional  values  of  g(XpXg)  except 
for  those  points  lying  on  the  line  X^  «*  3.  Portions  of  the  new  functions 
^(XpXg)  and  k(3^,Xg)  appear  in  Figure  Soane  functional  values  of 
k X^,Xg)  ore  in  square  brackets. 


Figure  b 


• ' ’ 


For  suitably  small  values  of  !,«/!  the  lattice  point  (X^Xg)  ~ (2,-1) 

has  a smaller  functional  value  than  any  of  its  adjacent  lattice  point.  , 


4 


vhich  by  the  definition  makes  (2P  -l)  a loe  1 minimum  for  the  integer- 

V 

valued  function  gtiL^Xg),  In  contrast,  the  continuous  function  fCx^Xg) 
does  not  have  a local  minimum,  only  an  absolute  minimum  at  (3,2}„ 

Regarding  round  off,  consider  any  point  on  the  diagonal  line  of 
in  Figure  2„  The  diagonal  line  is  a nonunique,  absolute  mini- 
mum for  f(X  ,X_)j  however,  for  many  points  along  the  diagonal,  rounding 

off  will  yield  an  integer  solution  which  is  clearly  not  the  solution 
to  the  integer  valued  problem  of  finding  the  minimum,  of  Only 

for  points  on  the  diagonal  with  1 > Xg,  or  Xg4  0 will  the  round  off 
yield  the  correct  answer „ 

We  have  demonstrated  that  even  with  a convex  function  having 
found  its  minimum,  rounding  off  to  an  adjacent  integer  point  does  not 
guarantee  a solution  to  the  corresponding  integer-valued  problem 0 
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